Circuit theory of unconventional superconductor junctions 
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We extend the circuit theory of superconductivity to cover transport and proximity effect in meso- 
scopic systems that contain unconventional superconductor junctions. The approach fully accounts 
for zero-energy Andreev bound states forming at the surface of unconventional superconductors. 
As a simple application, we investigate the transport properties of a diffusive normal metal in se- 
ries with a d-wave superconductor junction. We reveal the competition between the formation of 
Andreev bound states and proximity effect, that depends on the crystal orientation of the junction 
interface. 



In the last decade the mesoscopic superconducting 
systems have been the subject of intensive experimen- 
tal and theoretical research. The transport in these 
system is essentially contributed by the so-called An- 
dreev reflection a unique process specific for normal 
metal/superconductor interface. Each electron reflected 
from the interface may transfer a charge 2e to the super- 
conductor H, ||]. The phase coherence between incom- 
ing electrons and Andreev reflected holes persists in the 
normal metal at mesoscopic length scale. This results 
in strong interference effects on Andreev reflection rate. 
Q. The transport properties of mesoscopic N/S junc- 
tions have been theoretically investigated with various 
approaches, e.g., traditional nonequilibrium supercon- 
ductivity approach tunneling Hamiltonian approach 
Q , scattering formalism ||] and computer simulation . 

One of the authors has proposed a generic circuit the- 
ory of non-equilibrium superconductivity which accounts 
for the effects abovementioned. The mesoscopic system 
is presented as a network of nodes and connectors. A 
connector is characterized by a set of transmission coeffi- 
cients and can present anything from ballistic point con- 
tact to tunnel junction. Full isotropization of electrons 
is assumed in the nodes. This approach considerably 
simplifies a practical transport calculation, numerical as 
well as analytical. The circuit theory is based on conser- 
vation laws for so-called spectral currents. These addi- 
tional conservation laws present interference of electrons 
and holes. The spectral currents through each connector 
arc functions of spectral vectors in the nodes. There is 
one-to-one correspondence between spectral vectors and 
currents and Keldysh Green functions in the underlying 
microscopic approach. |^ Kirchoff-type equations deter- 
mine spectral currents and vectors in each node and con- 
nector, and, consequently, electric current in the circuit. 

Unconventional superconductors bring about very un- 
usual interface physics. The transport through the inter- 
face is influenced by formation of Andreev bound states 
(ABS) at this interface. 0, |ll|, |l2|. Those result 



from the interference of injected and reflected quasipar- 
ticles. The ABS manifest themself as a zero-bias peak in 
tunnehng conductance(ZBCP) [l2|. Indeed, ZBCP 



has been reported in various superconductors that have 
anisotropic pairing symmetry. [T^ The proper theory 
of transport in the presence of ABS has been formu- 
lated |]ll|, |l^ for conditions of ballistic transport only. 
This theory has to be revisited to account for diffusive 
transport in the normal metal. The point is that the 
diffusive transport provides an Andreev reflection mech- 
anism for ZBCP which does not involve any unconven- 
tional superconductivity. This mechanism may compete 
with the formation of ABS. The anomalous size depen- 
dence of transport in YBCO junctions reported in recent 
experiment seems to arise from this competition. 

All this has motivated us to extend the circuit theory 
to the systems containing unconventional superconductor 
junctions. We stress that this extension is by no means 
straightforward. The circuit theory can not be directly 
applied to an unconventional superconductor since it re- 
quires the isotropization. The latter is just incompatible 
with mere existence of unconventional superconductiv- 
ity. Fortunately, there is a way around. We concentrate 
on the matrix currents via the unconventional supercon- 
ductor junction to/from diffusive parts of the system. If 
one knows the relation between these currents and the 
spectral vectors (isotropic Green functions) in the diffu- 
sive part, one is able to use Kirchoff rules to complete 
the evaluation of the matrix currents everywhere in the 
system. 

This relation shall be derived from microscopic the- 
ory and presents the main result of this work. We stress 
that applicability of this relation is not restricted to cir- 
cuit theory. One can regard our result as a boundary 
condition for the traditional Keldysh-Usadel equations 
of non-equilibrium superconductivity. As an immedi- 
ate application, we study a d-wave superconductor junc- 
tion in series with normal metal. The resistance of the 
system appears to depend strongly on the angle a be- 
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tween the normal to the interface and the robe direction 
of d-wave superconductor (misorientation angle). This 
reveals the competition between the effect of ABS and 
proximity-induced reflectionless tunneling. 

To derive the relation between matrix current and 
Green functions, we make use of the method proposed 
in 1^. The method puts the older ideas |^ to the frame- 
work of Landauer-Biittiker scattering formalism. One 
expresses the matrix current in a constriction in terms of 
one-dimensional Green functions gn,tT;n' ,<t' (e; x, x'), where 
n,n' and a, a' = ±1 denote the indicies of transport chan- 
nels and the direction of motion along x axis, respec- 
tively. The " check" represents the Keldysh-Nambu struc- 
ture. These Green functions are to be expressed in terms 
of the transfer matrix that incorporates all information 
about the scattering, and asymptotic Green functions 
Gi^2 presenting boundary conditions deep in each side of 
the constriction. The isotropization assumption requires 
that these G do not depend on channel number. Under 
this assumption, the current is universal depending on 
transmission eigenvalues only. Although the isotropiza- 
tion assumption is good for conventional superconduc- 
tors and normal metals, it fails to grasp the physics of 
unconventional superconductor where the Green function 
essentially depends on the direction of motion and thus 
on channel number. To avoid this difficulty, we restrict 
the discussion to a conventional model of smooth inter- 
face, assuming momentum conservation in the plane of 
the interface. Within the model, the channel number 
eventually numbers possible values of this in-plane mo- 
mentum and the transfer matrix becomes block-diagonal 
in the channel index. We thus solve Green functions 
{£,x,x') separately for each channel. Asymp- 
totic Green function in the unconventional superconduc- 
tor does depend on the direction of motion cr, 

G2;n,a;n,a ^ G'^_^^— hG^J— ^ (1) 

reflecting different asymptotic conditions for incoming 
(G^"') and outgoing (G^"^) wave in each channel. The 
asymptotic Green function Gi in normal metal is the 
same for both waves and all channels (see Fig. 1). All 
these matrices satisfy unitary relation (Gji')^ = G^ = 1. 

After some algebra we obtain the matrix current in the 
following form 

I=^T.lG,,B,n], (2) 

rn 



H^^^ = {Gif ±G'i!l^)/2. 

Here Em = Jm/(1 + \/l ^ Tm)^ is related to the trans- 
mission coefficient r,„ in a given channel m. The above 
relation reduces to isotropic result of Ref. M provided 
G2^ = G2_ = G2. The above 4x4 matrix relation is 
the main result of the present work, ft incorporates the 
most general situation and allows for many applications 
that involve unconventional superconductors. Below we 
provide a simple but extensive application example that 
both illustrates circuit theory method and demonstrates 
an interesting interplay of ABS and proximity effect. 

The circuit is the one given in Fig. f : diffusive conduc- 
tor of resistance Rd in series with unconventional super- 
conductor junction. We disregard decoherence between 
electrons and holes in the diffusive conductor, ("leakage" 
current in terms of Ref.|^), this is justified at energies 
not exceeding Thouless energy of this piece of normal 
metal. We restrict our attention to d-wave supercon- 
ductor, being the most practical example of the singlet 
unconventional superconductor that preserves time re- 
versal symmetry. For simplicity, we have in mind a 
"two-dimensional" superconductor made from the layers 
stacked in z-direction. z-axis lies in the plane of the in- 
terface and is normal to the plane of Fig. 1. The interface 
normal (x-axis) makes an angle a with the main crystal 
axis. The propagation directions of the waves are thus 
in xy-plane and are parameterized by the angle 9 with 
X-axis. The angular dependence of the superconducting 
order parameter is thus given by A{6) — Aq cos(2(6'— a)). 
A scattering channel consists of an incoming wave in di- 
rection TT — 9 and outgoing wave in the direction 9. The 
sums over channels can be reduced to integrals over 9: 

2J oc / COS 61 (3) 

The Green functions are fixed in " US" terminal and in 
"N" terminal, the voltage V is applied to "N" terminal. 
The Green function in the node " DN" , Gi is not fixed and 
shall be determined from the balance of the matrix cur- 
rents. There is a natural separation of balance equations 
for 2 X 2 spectral currents that set advanced or retarded 
part of Gi and for particle current at a given energy that 
set the distribution function in the node " DN" . [|| 

We address the balance of the spectral currents first. 
The advanced 2x2 Green functions are fixed in "N" 
and "US" terminals and read Gat = r^, G2± = 

{A±Tx — itTz)/ ^ Aj_ — e^, T being Pauli matrices, A± — 
Ao cos(2(0±a)) being superconducting order parameters 
that correspond to direction of incoming (outgoing) wave. 
This suggests that the corresponding Green function in 
"DN" node assumes a form sin 7 • + cos 7 • where 
7 is yet to be determined. 7 is the measure of proxim- 
ity effect. All spectral currents are proportional to Ty. 
The spectral current ij^ through diffusive conductor is 
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proportional to the spectral angle drop,(|^), the spectral 
current i^-* via the interface is obtained from Eq.||. The 
balance equation thus reads 



-4^) = o 



(4) 



2e2 



E 



Under conditions considered, the transport is deter- 
mined by energy-symmetric distribution function, that is 
conventionally called /(. [|j The balance of particle cur- 
rents at each energy determines this distribution function 
in "DN" node. We will assume that the temperature \/ (3 
is much smaller than the typical value of superconducting 
energy gap, so we can disregard quasiparticle excitations 
in the superconductor. The particle current through dif- 
fusive conductor is given by the drop of the distribution 
function at its ends, the particle current via the interface 
is given by the corresponding block of Eq. ^. This yields 

it^ + 4^^ = (5) 

{ft - k)/RD 



;(P) 



2e 



(p) 

D 



fo being the symmetrized distribution function in the 
normal reservoir, fo = i[tanh(/3(e + ey)/2) — tanh(/3(e — 
eV)/2)]. The above relation becomes especially trans- 
parent if one regards T*'s as effective transmission co- 
efficients in each channel. It just shows that the full 
(energy-dependent) resistance of the system is the sum 
of the resistance of diffusive metal and the interface re- 
sistance, the latter being influenced by proximity effect. 
The degree of proximity effect is determined from Eq. 
(^. If we define the average over the angle as 

<A{e)>= ddcos9A{e)/ ddT{d) COS 9 

J-tt/2 J-tt/2 

with T{9) = Tjn, both balance equations can be rewritten 
in a compact form. 



R = Rd + RbI <T*{^,e,Tr„) >, 
7=<F(7,e,T™)>i?i3/i?B. 



(6) 
(7) 



Here i?s is the interface resistance in normal state, R is 
the full resistance. It may depend on energy, so that the 
full electric current is given by el^i — J defo{e)/R. 

To reveal the underlying physics, we present the 
concrete expressions for F = F{'j,e,Tjn) and T* = 
T*{'y,e,Tm) assuming | e |<c| A± |. It turns out that 
these expressions are essentially different for A+A_ >< 
0, this manifesting the formation of ABS in the latter 
case. For A+A_ < (ABS channels) we have: 



F 



-2T„ 



sm7 



Tm cos 7 - z(2 ^ 
T2(H-|cos7p 



- I sin7 |2) 



r2cos2 7+(2-T„)2(e/A)2 



-2tan7 (8) 



(9) 



with A = (2 I A+ II A_ |)/(| A+ I -h I A_ I). It is some- 
what counterintuitively that the zero-energy limit does 
not depend on the actual transmission, giving finite cur- 
rents even for insulating interfaces. This is the signature 
of the resonance forming precisely at zero energy [ pT| . 
If the transmission is low, the resonance feature persists 
in a narrow energy interval ~ TmA± only. The spec- 
tral current F eventually suppresses the proximity effect. 
The explanation is that ABS form a reservoir of normal 
electrons within the unconventional superconductor, and 
F can be viewed as a connection to this normal reservoir. 
The effective transmission coefficient T* at resonance is 
always bigger than 2, and is enhanced by proximity ef- 
fect. One can understand this as a mupliple Andreev 
refiection induced by the corresponding ABS. 

In the case of A_|_A_ > ("conventional" channels) 
the resonance feature is absent and energy dependence 
can be safely disregarded. The expressions are identical 
to those of conventional superconductor 



F = 



2TmSCOS7 



2 - T„i + TmS sin 7 
2r™[r™ + (2-r™)ssin7] 



2-T„ 



sin 7 



(10) 
(11) 



cos^ 7 



Here s = sgn(A+) = sgn(A_). The spectral current F 
thus induces proximity effect of the corresponding sign s. 
The effective transmission T* does not exceed 2 (which 
is the limiting case of ideal Andreev reflection). Being 
compared with the transmission in the normal state, the 
effective transmission is suppressed (enhanced) at Tm < 
(>)2/3. The fully developed proximity effect (7 — sn/2) 
restores the normal transmission. 

To summarize, the proximity effect originates from the 
"conventional" channels and is suppressed by ABS chan- 
nels. While the proximity effect is present, it enhances 
transmission via ABS channels. It restores the effective 
transmission of "conventional" channels to that in the 
normal state. The full resistance of the structure is de- 
termined by competition of all these effects. It is essential 
that one can tune the relative number of " conventional" 
and ABS channels by changing the misorientation angle 
a. As one can see from the Fig.l, the ABS channels are 
there in the angle interval tt/4 — \a\ < \6\ < n/A + \a\. 
If a = 0, there are no such channels. If a = 7r/4, there 
are no "conventional" channels. This gives no chance to 
proximity effect. 

To illustrate this further, we calculate with Eqs. (^,(0) 
the zero-voltage resistance (e 0) at different values 
of a as a function of Rb/Rd- The angular depen- 
dence of the transmission coefficient was assumed to be 
T{9) = cos^ 9/{cos^ 6+Z) with barrier parameter Z. The 
results are presented in Figs. 2,3. At Rd — there is no 
proximity effect in " DN" and the resistances are given by 
the quasi-ballistic formulas of Ref. jll]. The proximity 
effect may develop with increasing Ro and decreases the 
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interface resistance. This gets the curves down. The 
curves (a) and (a') correspond to d-wave junction at 
a — and conventional superconductor junction, respec- 
tively. One sees that the proximity effect is weaker in 
d-wave system. This is due to competition of the "con- 
ventional" channels having different signs of A+. The 
ABS channels appear with increasing a. The curve (6) 
demonstrates interface conductance reduced slightly be- 
low its normal state value. This manifests the enhanced 
transmission in ABS channels. The ABS channels quench 
the proximity effect very efficiently at a > 0.027r. The 
total conductance can be approximated hy Rd + Rro=o, 
although this becomes exact only at a = 7r/4. Follow- 
ing Ref. jl^, we regard the counterintuitive negative 
sign of [dR/ dRD)RD=o as a signal of importance of the 
proximity effect [or reflcctionless tunneling (RLT) ]. We 
evaluate this sign at different Z and a. (Fig. 3) The sign 
of dR/dRo is negative for junctions of low transmissivity 
in a relatively narrow range of a. 

In conclusion, we have extended the circuit theory of 
superconductivity to include unconventional supercon- 
ductor junctions. We have derived a general relation 
for matrix current to/from unconventional superconduc- 
tor. An elaborated example demonstrates the interplay 
of ABS and proximity effect in a c?-wave junction. The 
theory presented will facilitate the analysis of more com- 
plicated mesoscopic systems that include unconventional 
superconductors. 

The authors appreciate useful and fruitful discussions 
with J. Inoue, H. Itoh and G.W.E. Bauer. 




FIG. 1: The unconventional superconductor junctions (black 
box) can be incorporated into circuit theory by means of 
the matrix current relation (^). This relation accounts for 
anisotropic features of the US, as sketched for a d-wave su- 
perconductor. 
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FIG. 3: "RLT" ("NRLT") marks the region where 
dR/dRo < {dR/dRo > 0) at Rd = 0. Tr is the 
average transmissivity of the junction in the normal state. 
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